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[57] ABSTRACT 

Topographical terrain models are generated by digitally 
delineating the boundary of the region under investiga- 
tion from data obtained from an airborne synthetic aper- 
ture radar image and surface elevation data concur- 
rently acquired either from an airborne instrument or at 
ground level. A set of coregistered boundary maps thus 
generated are then digitally combined in three dimen- 
sinoal space with the acquired surface elevation data by 
means of image processing software stored in a digital 
computer. The method is particularly applicable for 
generating terrain models of flooded regions covered 
entirely or in part by foliage. 

17 Claims, 8 Drawing Sheets 
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GENERATION OF TOPOGRAPHIC TERRAIN ouiviMARY 

MODELS UTILIZING SYNTHETIC APERTURE It is an object of the present invention, therefore, to 

RADAR AND SURFACE LEVEL DATA combine synthetic aperture radar data with other 
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ORIGIN OF THE INVENTION 

The invention described herein was made by an em- 
ployee of the United States Government and may be 
manufactured and used by or for the Government for 
governmental purposes without the payment of any 
royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

This invention relates generally to terrain mapping 
and more particularly to the generation of multi-dimen- 
sional topographical terrain models of the earth’s sur- 
face such as forested regions which have an underlying 
surface which is periodically flooded or includes a 
water surface of varying elevation. 20 

Wetland forests and Mangrove or wet coastal forest 
ecosystems are an important resource throughout the 
tropic and subtropic regions of the globe. As a whole, 
these formations account for a significant share of the 
forested coast lines and river flood plains of these re- 25 
gions and are an important source of biological produc- 
tivity. By providing environments for abundant plant, 
algal, and marine and terrestrial animal life these ecosys- 
tems figure Prominently in the cycling of organic mat- 
ter. 30 

A critical environmental factor in the formation, 
health, and status of these ecosystems is the surface 
hydrology, the two most important factors of which are 
the water level dynamics and water chemistry. Both of 
these factors are critically interactive with the geomor- 35 
phology or substrate topography upon which flood 
plain forests and mangroves evolve. The topographic 
contour of the floor of these forests, moreover, deter- 
mines the frequency of the flooding and the chemistry 
of the flood waters and thus the productivity and health 
of the forest and the utility of the soils. 

Heretofore, generation of contour maps or obtaining 
terrain information of forested regions relied solely 
upon leveling survey measurements acquired on the 45 
ground. This meant that a team of surveyors actually 
had to enter the area and take required measurements. 
Aerial survey was impossible as the terrain is hidden by 
the forest canopy. Ground survey is a very tedious and 
expensive process that requires total and long term 
access to the area of interest. Also, many of the world’s 
coastal forests, river deltas and flood plains are not 
readily accessible. These include the mouths of tfie 
Ganges river, the Amazon basin, the Mississippi delta, 
and the Mekong delta, to name a few. In addition to 55 
being very time consuming, such efforts are extremely 
tedious, require many people, extensive travel and logis- 
tics, and often access to the areas to be surveyed is 
impossible because they lie in a politically sensitive area. 

While terrain modeling and foliage penetration using 60 
a synthetic aperture radar is generally known and can 
also be used to remotely determine water level given 
the use of special “in situ” devices, radar derived flood 
boundary maps combined with tide and hydrology data 
is not known to have been used to generate topographic 65 
models of subcanopy surface areas which periodically 
change due to the water surface variation caused bY the 
changes in the water level in the region of concern. 


known surface information to generate a topographic 
image of a geographical surface area. 

It is a further object of the invention to combine 
synthetic aperture radar images of the earth’s surface 
obtained from an airborne radar with concurrently 
acquired tide surface information to generate a digital 
terrain model of a predetermined water surface area. 

It is another object of the subject invention to utilize 
synthetic radar imagery for generating topographic 
information of a subcanopy flooded surface area by 
combining the radar derived boundaries of inundated 
water surface areas with concurrently obtained surface 
elevation data. 

And yet another object of the invention is to generate 
a multi-dimensional terrain model of a certain region of 
the earth for use in determining tidal flow and eventual 
nutrient transport from a forested area to a marine or 
river ecosystem by combining data acquired by syn- 
thetic aperture radar making periodic passes over the 
area under investigation and combining it with known 
water levels acquired either from the ground or from 
above. This can also be used to monitor flood damage 
and the effect of ground subsidence or rising sea levels. 

Briefly, the foregoing and other objects are achieved 
by a method and apparatus for generating topographi- 
cal terrain models by delineating the boundary of the 
region under investigation by data obtained from an 
airborne synthetic aperture radar. The radar data from 
an airborne radar platform is thereafter combined with 
concurrently obtained surface information to generate a 
digital terrain model via a digital computer. Of particu- 
lar concern are forested regions having an underlying 
variable water surface area, but of equal importance are 
river flood plains of low relief and other types of vege- 
tated areas. The dominance of the interaction phenome- 
non, the canopy to surface or surface to canopy reflec- 
tion in a flooded forest is utilized to create subcanopy 
flood boundary maps for at least two different time 
periods. The boundary maps are then digitally com- 
bined in three dimensional space with acquired water 
surface elevation data by means of image processing 
software stored in a digital computer. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The following detailed description of the invention 
will be more readily understood when considered in 
conjunction with the accompanying drawings wherein: 

FIG. 1 is a diagram illustrative of the manner in 
which radar data is acquired in accordance with the 
subject invention; 

FIG. 2 is a diagram illustrative of radar flood bound- 
aries acquired at three different time intervals for a 
predetermined region; 

FIG. 3 is a diagram illustrative of the merging of the 
radar derived flood boundaries shown in FIG. 2 and 
concurrently acquired water elevational data for gener- 
ating topographic contour map in accordance with the 
subject invention; 

FIG. 4 is a diagram further illustrative of the merging 
of radar derived flood boundaries with water eleva- 
tional data;. 

FIGS. 5A and 5B are diagrams illustrative of the 
preferred method of the subject invenifon; 
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FIG. 6 is a diagram further illustrative of a portion of 
the method shown in FIGS. 5A and 5B; and 

FIGS. 7A and 7B are illustrative of a flow chart of 
computer software for implementing the method shown 
in FIGS. 5A and 5B. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Referring now to the drawings wherein like refer- 
ence characters refer to like elements, FIG. 1 discloses 
a radar platform 10 which may be, for example, a space 
shuttle, airplane, or satellite which carries a synthetic 
aperture radar (SAR), not shown. Synthetic aperture 
radar comprises apparatus well known to those skilled 
in the art and preferably comprises an L-band, P-band 
or C-band radar having quad polarization or any combi- 
nation of HH, VV, HV, VH, RC, LC polarizations and 
which is capable of penetrating vegetation, particularly 
tree foliage and which operates with an incidence angle 
which ranges between 25° and 57° . The radar generates 
an image, preferably but not necessarily, at the same 
incidence angle of the earth’s surface along a relatively 
narrow strip 12 during a pass over the area of interest 
during an orbital flight. As shown in FIG. 1, the flight 
path proceeds from the Bay of Bengal to the mouth of 
the Ganges river where it passes over Bangladesh and 
then continues northward over China. Reference nu- 
meral 14 comprises a survey site located, for example, at 
the mouth of the Ganges river in southern Bangladesh. 

In the subject invention, SAR surface image data is 
acquired over the survey site 14, typically but not re- 
stricted to, 70km by 40 km, at predetermined time inter- 
vals, for example, different stages of inundation of a 
forest region caused by flooding. The period of time 
involved, however, is selective and can be hours, for 
example, for tides and days or months and even years 
for seasonally flooding rivers. Simultaneously with the 
SAR data acquisition, ground water level measure- 
ments are acquired in situ from a plurality of gauges or 
other type of depth measuring devices located at one or 
more, typically at least five stations 16 within the survey 
area, or it is acquired remotely by means of radar altim- 
etry. When desired, conventionally acquired leveling 
survey transect data within the survey area and ob- 
tained by making elevation measurements every 30 
meters, for example, for a distance of 1200 meters can 
also be incorporated into the model generation in order 
to obtain enhanced accuracy. The water elevation data 
thus acquired are then entered into a digital computer 
where stored routines, hereinafter described, generate 
one or more surface models. The radar images are then 
geometrically matched to the water surface models by 
being digitally coregistered therewith using, typically, a 
first order polynomial interplation algorithm by means 
of the computer software from which three dimensional 
surface models are generated. 

The method of the subject invention includes three 
major processes which are: digital filtering, bit masking 
and interpolation. The radar data is converted to statis- 
tical data for a plurality of test sites on the ground 
where it is processed digitally through the use of spatial 
filters to generate one or more relative brightness maps. 
The maps are divided into brightness classes with the 
division being made based on either a statistically signif- 
icant break point or on ground data indicating a place 
which corresponds to a place on the SAR image where 
an actual flooded or non-flooded boundary is located 
within the survey area. This division then results in a 
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flood boundary map being generated for the area under 
consideration and is performed for each radar image 
generated. The flood boundaries for each image data set 
are next stripped off using bit masking or a boundary 
5 tracing algorithm. These boundaries are then multiplied 
to the appropriate water elevation surface models to 
generate real elevations from the boundaries. The new 
elevation data from one or more such boundaries are 
then compiled into a single composite image map. An 
10 interpolation algorithm is then performed on the com- 
posite elevation data to generate a digital three dimen- 
sional elevational model. With this technique, contour 
maps, perspective views, or a variety of volume and 
topographic products are selectively generated. 

15 This procedure is illustrated in FIGS. 2, 3 and 4. 
Referring now to FIG. 2, shown thereat are three flood 
boundaries generated from SAR radar image data of the 
survey site 14, for example at three different time inter- 
vals Ti, T2and T3. At the time Ti, a flood boundary A 
20 is produced, whereas at the time T2, a flood boundary B 
is generated which is relatively smaller than the flood 
boundary A. This is indicative of less flooding in the 
area at the time T2. At the time T3 f the flood waters 
have receded even more, generating a flood boundary 
25 C which is shown in relation to the largest flood bound- 
ary A. 

In the invention, the three flood boundaries A, B and 
C generated are mutually centered as shown in FIG. 3 
on the centroid D where they are merged with water 
30 surface models A', B' and C' generated from the ground 
surface data concurrently obtained at a plurality of test 
sites or stations 16 at the times Ti, T2 and T3 which 
coincide with the times the radar images for the flood 
boundaries A, B and C were generated. The merged 
35 image data is then used to generate a three dimensional 
topographic terrain model of the region 14 or a portion 
thereof as shown, for example, in FIG. 4. 

Referring now briefly to FIG. 4, shown thereat are 
two radar images 18 1 and I82 which are generated at the 
40 times Ti and T2, respectively. A plurality of test sites 20 
are located on the SAR radar images 18 1 and I82 from 
which statistics are generated. The radar images are 
then spatially filtered by median or other statistics based 
filters 22 1. and 222. Relative brightness maps 24 1 and 242 
45 are next generated for the times T 1 and T2. The brigness 
maps 24i and 242 are partitioned into brightness classes 
at 26 1 and 262 in response to the statistics acquired from 
the test sites 20 from which the flood boundary images 
A and B (FIG. 2) are generated. These boundaries are 
50 then merged with the water surface models A' and B' to 
generate a final composite three dimensional terrain 
model 28. 

Further details of the method and the computer soft- 
ware for implementing the topographic terrain model 
55 generation will now be considered. As depicted in 
FIGS. 5A and 5B, the initial steps involve mutually 
separate processes as they relate to the SAR image data 
and to the water elevation data. At three different times 
Til, T2 and T3 } the SAR acquires three separate data sets 
60 over the survey area 14 as indicated at 30A. The SAR 
operates in a frequency range which has vegetation 
penetration properties. The data is collected in or con- 
verted to the digital domain for image manipulation 
purposes. Although one set of data for a time T n will 
65 suffice in certain instances, two or more data sets are 
preferred and are acquired at times where different 
water levels due, for example, to flooding, tides or sea- 
sonal hydrological regimes are in effect. As noted ear- 
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Iier, the time interval between SAR data acquisitions 
can be hours, months, or even years. The radar images 
are also geometrically coregistered so that a statistical 
analysis can be performed on a single geographic defini- 
tion of a test site in the areas on the imagery. 

The initial step 30A of acquiring the SAR data is 
followed by two steps 32A and 34A where flood bound- 
ary delineation and flood boundary images are gener- 
ated. These two steps are shown in further detail in 
FIG. 6. 

Referring now briefly to FIG. 6, step 32 A i designates 
a preliminary step used for providing accurate flood 
boundary delineation and comprises what is referred to 
as the acquisition of ground truth or on-site observa- 
tions. Although not always necessary, ground truthing 
provides a higher level of confidence in the accuracy of 
flood boundary delineation on the SAR images gener- 
ated. To implement ground truth, a predetermined 
number of test sites 20 are located on the SAR image 21 
where sub-canopy or flood plain conditions can be veri- 
fied on the ground, i.e., flooded vs. non-flooded states 
during SAR data acquisition. In forested areas, the test 
sites 20 are selected that represent the major forest types 
and in topographic areas that will yield flooded and 
non-flooded examples for each set of SAR data. Statis- 
tics are then generated for digital analysis of the sites 20 
in step 32A2. Following step 32A2, several techniques 
can be utilized to define flood boundaries, one being 
controlled density slicing. Another is the use of polar- 
ization signature analysis to determine the presence of a 
two or multiple bounce signature (flooded area) verses 
a single bounce (non-flooded area) signature. This ob- 
tains for a flooded canopy; however, these relationships 
may be reversed for non-forested areas and open water. 

Controlled density slicing involves a statistical analy- 
sis of the SAR image pixel values within the designated 
ground truth test sites 20 and includes means for deter- 
mining: means, medians, standard deviations, etc. Fol- 
lowing this, the data is subjected to spatial filtering as 
shown in step 32Aa of FIG. 6. Since vegetated areas are 
of concern, a median or mean guided spatial filter is 
provided to fit a predetermined regional scale such as 
5x5, 7x7, or 11x11 meters, where each element is 
approximately 25 meters square for forested areas with 
relatively homogeneous stands and applied to the SAR 
data sets. In any event, the size of these filters are deter- 
mined by the size of the forest and the resolution of the 
radar. This yields a fine density pattern over the image. 
Statistics can again be regenerated on the sites 20 after 
filtering for a simpler approach to determining flood 
boundaries. 

Next a boundary point selection and gradient remap- 
ping step 32 A4 is performed using a gray scale separa- 
tion using the test site statistics generated either in step 
32A2 or 32A3. The flood boundary is delineated as being 
the mid-Point value between flooded and non-flooded 
grouped means or medians. A generic gray scale density 
slicing routine is used to separate classes into flooded vs. 
non-flooded categories by remapping the images into 
two classes above and below a statistical boundary 
Following this, a boundary isolation step 34Ai is per- 
formed where the boundaries are prepared for subse- 
quent merger by bit masking, i.e. the flood boundary 
image pixels are assigned the digital values of 1 and all 
others 0. An error checking step 34A2 is also shown. 
This is performed where any values greater than 1 indi- 
cate overlap Problems and the step 32A4 is repeated, if 
necessary. Following the boundary isolation and bit 
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masking step of 34A, the flood boundary images are 
merged with water elevation values as shown in step 36 
of FIG. 5. 

Referring now back to FIG. 5A, step 30B is intended 
5 to designate that water elevation data is also being ac- 
quired over the survey area at the same time that SAR 
data is being obtained as shown in FIG. 1. The location 
of the water level measurements with respect to the 
survey site and the SAR image data is known and com- 
10 prises the sites 16 in the survey area 14. In order to 
generate a first order trend surface, three or more sites 
distributed over the survey site are preferred. The 
water level data is acquired and extracted to match the 
time Tj, T2 and T3, for example, of the SAR image data 
15 acquisition shown in step 30A. 

Next step 32B is performed, which comprises the 
generation of water surface models. To accomplish this, 
water elevation data is converted into a digital image 
format to form a grid. Elevation values are then interpo- 
20 lated by an interpolation routine, with elevation images 
being generated as shown by step 34B, which is illustra- 
tive of the elevation of the water surface over the entire 
area. This procedure is applied to each elevation data 
set corresponding to the SAR data acquisition sets and 
25 are registered geometrically with the SAR generated 
images relative to a common base, the common base 
being, for example, the centroid D shown in FIG. 3. 
Any standard or custom map projection, however, can 
be utilized, when desired. 

30 With the generation of the water surface elevation 
models, a merger of the SAR flood boundaries with the 
water surface elevation models is shown by step 36, 
whereby each SAR boundary set is converted to real 
elevation values by matching it to its appropriate water 
35 surface elevation model. This is accomplished by 
known digital data processing techniques. For example, 
the results obtained by the flood boundary generation 
step 34A is multiplied in the image domain with the 
respective counterparts of the water elevation models 
40 generated in step 34B. The SAR flood boundaries 
which have a digital value of 1 will acquire the appro- 
priate elevation values from the water elevation image. 
All other values being 0 will retain the value of 0. Step 
36, moreover, is repeated for all SAR flood boundaries 
45 and water surface images over the desired time inter- 
vals, for example, the times T 1, T2 and T3. It should also 
be noted that the interface between land and water for 
islands can also be acquired in the same way. 

Referring now to FIG. 5B which comprises an exten- 
50 sion of FIG. 5A, reference numeral 38 is indicative of 
the step wherein the results of steps 36 and 37 are com- 
bined into a composite data set from which a single 
composite image 40 can be generated. The new data set 
comprises image data containing real elevation points 
55 over the survey area based on the SAR flood boundary 
and elevation data set merger. Error checking can fur- 
ther be utilized, when desired, to eliminate bad data. An 
interpolation algorithm is further performed on the 
results obtained in steps 36 and 37 in the generation of 
60 the composite image as shown by step 42. The result 
obtained is a complete terrain image of the survey site 
based on the SAR and water elevation data acquired 
over the time intervals selected. 

Reference numeral 44 indicates a step of product or 
65 image generation which can take many forms, several 
of which are shown comprising contour maps 46, per- 
spective views 48, elevation class maps 50 and inunda- 
tion frequency maps 52, to name a few. Such maps can 
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be used for forest productivity modeling in forest areas, 
development planning, disaster management and hydro- 
logical modeling, etc. 

Since any type of SAR polarization, e.g. HH, HV, 
VV or VH can be used, when desirable, polarization 5 
signature analysis may be used for defining the flood 
boundaries. This technique involves the identification 
of flood boundaries by verification of a single bounce 
vs. a two-bounce polarization signature. All zones ex- 
hibiting a strong two bounce return can be considered 10 
flooded in vegetated areas. A specular or single bounce 
may indicate open water in non-vegetated areas. Once 
the signature has been determined, these areas can be 
classified by well known generic classification proce- 
dures. The flood boundary is then determined to be an 15 
image element where a strong two bounce signature 
gives way to other non-indicative signatures, i.e., any- 
thing other than two bounce signatures. 

The method outlined in FIGS. 5A and 5B as well as 
FIG. 6, can, when desirable, be implemented by a soft- 20 
ware package entitled, “Land Analysis System (LAS) 1 ’ 
developed by and for the Goddard Space Flight Center, 
National Aeronautic and Space Administration. This 
computer software includes a series of programs or 
routines which digitally operate on both the SAR data 25 
and the water elevation data acquired in order to gener- 
ate the required terrain models from the water surface 
elevation models and SAR derived flood boundaries. 

Referring now to FIG. 7A, reference numeral 54 
designates five programs which operate on the acquired 30 
SAR data for coregistering the SAR data sets to a com- 
mon geographical base. These programs are shown 
entitled: TEIPTS, TEIMERGE, NULLCORR, TIE- 
FIT and GEOM. Next the data is fed to a STATS 
program 56 to delineate the test sites on the SAR images 35 
generated. Next a MIN-MAX program 58 generates 
statistics for the test sites delineated. This is followed by 
a two program set 60 identified as LOWCAL- 
MEDIAN and LOWCAL-MEAN which spatially fil- 
ters the radar data This is followed by boundary selec- 40 
tion using the results from the MIN-MAX program 58, 
which can be reapplied to the process at 62, if desired, 
and the boundary gradients are defined by a RENUM- 
BER routine 64. The boundary locations on the images 
are next generated using the program identified as 45 
UGRAD 66 which is again followed by another RE- 
NUMBER routine 68 which bit masks the boundary to 
digital number Ts and 0’s. This data is now fed to a 
merging program 70 shown identified as MULTIPLY 
in FIG. 7B. An ADDPIC program 72 is also included as 50 
shown in FIG. 7A in order to perform error checking 
on the water boundary data generated at 68. 

Referring now to FIG. 7B, the water elevation data 
acquired is operated upon by two programs 74 identi- 
fied as TRANCORD and SURFACE-POLY. The 55 
TRANCORD program converts water level data to an 
image format registerd to the SAR generated images. 
The SURFACE-POLY program generates a water 
elevation surface model which is applied to the MUL- 
TIPLY program 70 where each SAR boundary image 60 
is merged with its appropriate water elevation model. 
Next the ADD-PIC program 76 is again run. This 
merges all of the results of the MULTIPLY program 70 
into one image set of real elevation values. Next an 
interpolation algorithm is provided by two programs 65 
identified by reference numeral 78 as ZOOM RE- 
DUCE and SURFACE-OCTANT, the first of which 
operates to reduce the data set to a manageable size, if 
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desired, while the second performs the interpolation 
and generates a digital elevation model of the survey 
set. 

Finally, a display of the elevation model in a relief 
map format or perspective view, is performed by a pair 
of routines shown by reference numeral 80 and identi- 
fied as MAP and PER SPEC. 

Thus what is shown and described is an extremely 
powerful tool which can utilize off the shelf image 
processing software to generate three dimensional ter- 
rain models which would be extremely difficult, if not 
impossible, to obtain with conventional apparatus and 
methods. 

Having thus shown and described what is at present 
considered to be the preferred embodiment of the in- 
vention, it should be noted that the same has been made 
by way of illustration and not limitation. Accordingly, 
all alterations, changes and modifications coming 
within the spirit and scope of the invention as set forth 
in the appended claims are herein meant to be included. 

I claim: 

1. A method of generating topographical terrain 
model of a flooded region of the earth’s surface cov- 
ered, at least in part, by foliage, comprising the steps of: 

(a) obtaining radar generated surface image data of 
said region; 

(b) generating a boundary map delineating two differ- 
ent surface types in said region for at least one time 
period from the radar generated surface data; 

(c) concurrently acquiring water surface elevational 
data of said region during said time period; 

(d) combining said boundary map with said water 
surface elevational data; 

(e) generating a multi-dimensional image of said re- 
gion from the combining of said boundary map 
with said water surface elevational data; 

said step (a) of obtaining radar generated surface 
image data and step (c) of acquiring water surface 
elevational data comprises obtaining and acquiring 
said data for a plurality of mutually concurrent 
time periods and wherein step (b) of generating a 
boundary map comprises generating at least one 
flood boundary map for each of said time periods; 

(0 geometrically coregistering the flood boundary 
maps generated in step (b) and then combining the 
coregistered boundary maps with said water sur- 
face elevational data acquired during said time 
periods; and, 

wherein said step (e) comprises generating a three 
dimensional water surface model of said region 
from said water surface elevational data. 

2. The method of claim 1 wherein said step (a) of 
obtaining radar generated surface image data comprises 
generating said data from an above ground means that 
could include airborne and spaceborne platforms. 

3. The method of claim 2 and wherein said step (c) of 
acquiring water surface elevational data comprises ac- 
quiring said data at ground level. 

4. The method of claim 3 and wherein said step (c) of 
acquiring water surface elevational data comprises ob- 
taining water elevational data at a plurality of sites 
within said region. 

5. The method of claim 2 and wherein said step (c) of 
acquiring water surface elevational data comprises ob- 
taining said elevational data by means which include 
ground level, airborne and spaceborne devices, and at a 
plurality of measurement sites within said region. 
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6. The method of claim 5 wherein said step (a) of 
obtaining radar generated surface data further com- 
prises the step of generating at least one synthetic aper- 
ture radar image of said region during each time period 
and converting said image into digital image data, 

wherein said step (c) of acquiring water surface ele- 
vational data of said one surface type further com- 
prises converting said data into digital data and 
generating a digital elevational water surface 
model therefrom, and 

wherein said step (d) of combining comprises digi- 
tally merging said radar image data with said eleva- 
tional water surface model. 

7. The method of claim 6 wherein said converting 
step of step (a) further comprises the steps of, 

(g) selecting a predetermined number of surface test 
sites on said radar image; 

(h) converting the image data of said test sites into 
statistical data; 

(i) spatially filtering said radar image data; and 

<j) generating at least one relative brightness map 
from the statistical data of said test sites from which 
said flood boundary map in step (b) is generated. 

8 . The method of claim 7 wherein step (b) includes 
digital bit masking of the image data and wherein sec- 
ond type surface data is strapped off for delineating the 
boundary of said two surface types. 

9. The method of claim 6 wherein said step (c) of 
generating a digital elevational water surface model 
includes performing a digital interpolation algorithm to 
generate a multi-dimensional elevational water surface 
model. 

10. The method of claim 6 wherein step (e) of gener^ 
ating said multi-dimensional image of said region in- 
cludes performing a digital interpolation algorithm on 
the merged radar image data and the water surface 
elevational data. 

11. The method of claim 7 wherein said one surface 
type comprises a water surface that can include a 
flooded land surface and the other surface type com- 
prises a relatively dry land surface that can include a 
non-flooded land surface, and 

wherein said boundary maps indicate respective land 
edge which includes flooded/non-flooded bound- 
aries therebetween 
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12. A method of generating a topographical terrain 
model of a flooded region of the earth’s surface cov- 
ered, at least in part, by foliage, comprising the steps of: 

(a) obtaining radar generated surface image data of 

5 said region; 

(b) delineating and generating data corresponding to 
the flood boundary of said region for at least one 
time period from the radar generated surface data; 

(c) concurrently acquiring water surface elevational 
10 data of said region during said time period; 

(d) merging the flood boundary data with the water 
surface elevational data; 

(e) generating a multi-dimensional terrain model of 
said region from the merging of the flood boundary 

15 data with the water surface elevational data; 

said steps (a) and (c) comprise obtaining and acquir- 
ing data for a plurality of concurrent time periods; 
step (b) includes the step: 

(0 of coregistering the respective flood boundaries 
20 for said time periods; and, 
additionally including the step: 

(g) of generating a three dimensional water surface 
model from the water surface elevational data ac- 
quired for said same time periods. 

25 13. The method of claim 12 wherein said step (a) of 

obtaining radar generated surface image data comprises 
generating said image data with a synthetic aperture 
radar. 

14. The method of claim 13 wherein said radar has a 
30 foliage penetrating operating characteristic. 

15. The method of claim 12 wherein said step (c) of 
acquiring the water surface elevational data comprises 
obtaining said data at ground level. 

16. The method of claim 12 wherein said step (c) of 
35 acquiring the water surface elevational data comprises 

acquiring said data from an above ground means that 
could include airborne and spaceborne platforms. 

17. The method of claim 12 wherein said step (a) 
further comprises the step of, 

40 (h) selecting a predetermined number of ground truth 

test sites within said region on said radar image; 

(i) converting the image data of said test sites into 
statistical data; 

(j) spatially filtering said radar image data; and 

45 (k) generating at least one relative brightness map 

from the statistical data of said test sites from which 
said boundary map in step (b) is generated. 

* * * * * 
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